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Summary 
The starting point is a function due to Dirichlet, the function f(s) 


—w Cb 
B(s) = Li n+ 


7 
If s = 7 the value is os by isolating 61 we find first of all an approximation of 61, indeed 


184320 
1= s— = 61.0271871... 
1 
Or more precisely 
256! 1 1 1 


So, by inverting this series we get a good approximation of a prime number and if we push further 
an equality. Here the number 61 is the 6 Euler number. The largest known Euler number that is 
prime is E;,) a number of 1062 digits, there does not seem to be any other before Ej 99999. We 


take here the numbers generated by the series ar where all the E,, are positive. It should also 


be noted that the approximation provided by this process gives us 


2512 510! 
E510 ~ qy511 


which is excellent since 243 exact decimals out of 1162. The goal of this article is to find all 
possible expressions that can represent a prime number. 


The Euler, Bernoulli and Tangent numbers 


As indicated in the summary, the Taylor expansion of the function 1/ 
cos (x) function provides the Euler numbers and if some of these are prime then we will 
have a representation of this prime. It remains to find which other primes can be 
represented. The Euler-Tangent numbers for example are obtained with the expression 


1 + sin (x) 
cos (x) 


Which once developed in exponential series generates the numbers 


1, 1, 1, 2, 5, 16, 61, 272, 1385, 7936, 50521, 353792, ... (sequence A000111 of the OEIS 


gnt2 ! 


“. We notice that 2 is a factor in the initial 


catalog) whose general term, a(n) ~ aa" 


. 272 32*7! 7 
expression, we can then represent 17 = qo ae well as 277 since a(8) = 1385 = 
1024+8! 8257536 , é . : . 
and thus 277 = We can then obtain an approximation of this prime as well 


™? nm 
as the exact expression by taking the series. Our list of primes is getting longer. 


The list so far is 2,5, 17, 31, 61, 277, 691 and E540. 


The expressions are always of the form 


c™k! 


p at (1) 


Where c and d are constants, in the examples above, c = 2 and d can be integer or 
algebraic. These first few are just the tip of the iceberg for an obvious reason. It is easy to 
generate numbers like Euler or Bernoulli with other trigonometric expressions like 


sin(x) + cos(3x) 
cos (4x) 
The development in Taylor series gives the following sequence: 1, 1, 7, 47, 497, 6241, 
95767, 1704527, 34741217, ... 
We see 7, 47. So if we can calculate the general term a(n) of this sequence we can add the 
first ones found to the list. Indeed, the general term is 
8"T(n) 


2V4+2V2 0" 


a(n) = 


That we find in this way. 


1) We develop in series the trigonometric expression 
2) We collect the coefficients with n! 
3) We calculate a(n + 1)/a(n). 

4) We calculate the differences term by term to find the constant 8/7. 
5) We go backwards (retro-engineering) to isolate the constant d in the expression 


(1). 


This makes it possible to obtain that 


83 


7 = ———., 47 


V4t2y2 73 


834! 


V4+2V2 14 


~ 
~ 


The approximations found are coarse when n is small but the 200" term is of the order 
of 10454 and a(200) has a precision of 95 decimal digits. 


If we look at other functions like the ones above, we find other prime numbers, 
here are some examples. All expressions are developed in exponential generating series. 
In bold, a term over 2 is considered. Annnnnn refers to the Online Encyclopedia of Integer 


Sequences: OEIS.org. 


Function Suite General term First Example 
A000464 
sin(x) 1, 11, 361, 24611, Ce Peace eee ee tw 2682 
cos (2x) 2873041, 512343611, = one ak oar 
1 —2cos (x) A000708 ae 
aE? -1, -1, 0,1, 6, 29, 150, 841, 7 3, 29, 
Ten) 5166, 34649, m 
cos (x) A000281 V2 2” (2n)! 42* 3 19 ax Siz 
cos (2x) 1, 3, 57, 2763, 250737, = oo ~ 78 
1 A000364 2a 28 6! 
cos (x) 1, 5, 61, 1385, per 561,277, | 61x = 
sin(x) + cos(3x) A006873 8"l(@) ye _ 8 
cos (4x) 1,1, 7,47, 497, 6241, 2J4+2V2 7" toy © 4422 03 
1+sin (x nt2 yI 4 
(x) A000111 oe al 25,17, 7 w 327! 
cos (x) 1,1, 1,2, 5, 16, 61, 272... qun+t 18 
cos(x) + cos(2x) + cos (3x 2n —?)! 52571 = 
EOS SOS) 3,11, 29, 191, 871, 52571, BEE a) spe wiees 212 10! 
cos (x) 7 (2n-1) 191, 52571 — 


cos(2x) + sin(2x 4n _4)1 
(2x) (2x) A012393 24" (2n — 1)! 17,31,691 Beas ceed 
cos (2x) 2,16,512,34816,4063232, et 72 


We could continue these examples to find other prime numbers, these few functions are 
among the simplest. In addition to isolated primes if a term is composed of 2 factors, it is 
always possible to obtain a new identity as in the case of the 8" Euler number : 1385 
which is 5 - 277. Since we already have 5, this provides an approximation for 277. The 
process works but it is a bit tedious to go and get any prime and illusory to find an 
expression for each prime. For the moment, this short list is only a sequence of exotic 
expressions giving almost a prime number. 


Other forms giving prime numbers 
Ramanujan found this formula with Bernoulli numbers and an unusual sum. 
AK+1 


n Bante 
» eo Ss (2) 
e 1 2(4n+2) 


n=1 


This gives for 4k + 1 = 13the identity 


If k is much higher, we have 
a no73 
24.) sax 7 = 1563446 ... 036059151 
n=1 


A prime number of 1077 digits, with the exponent 22607 we have one of 71299 digits, it 
is the largest known of its kind. What is remarkable for this sum is the affinity with once 
again the number zz. Remarkable also because of the exceptional precision. 


Indeed, it is not difficult to establish that in general we will have this 
approximation 


= k! 
>» erm _ 1 oe (rm)k*1 (3) 


n=1 


So, if the result of the infinite sum is prime we will have another representation of these 
approximations. If kis high, the approximation is spectacular, here with the exponent 673 
has 202 exact decimals. 


These sums (2) are known, what is less known is to extend the Ramanujan formula to 
two terms. 


Using an example with 691 we have (691 is the numerator of the 12‘ Bernoulli 
number). 


i. lt = nit 
691=16) sa —q-2) samy (4) 
n=1 n=1 
The first sum already provides a good approximation of 691 2 Te = - Again, this 


approximation can be made exact by adding either the missing term or a as of the series 
of the 2"4 member of the equation. This equation has a strong resemblance to one given 
by Jean-Pierre Serre (A course in arithmetic p.157). 


65520 g,4(n)q” 

11(1)q 

691 (>) 
n=1 


In disguise the series (2) are Eisenstein series where the variable q” becomes e?”". The 
identity (4) with 691 is just a rewriting of (5). 


Eisg = 1+——— 


A summary analysis of the identity (2) and (3) shows that we have either 1 term or 2 
terms which give a rational, for example 


1=24) 4-264) aang (6) 
n=1 n=1 


We can obtain an infinity of them as (6) if the exponent is of the form 4k —1. The 
coefficients grow with k, when k = 5 we have 


9 19 
aa > e4mn _ 1 
n=1 n=1 
(7) 
1946319750384384375 
221930581 ~= ==. = 221930369.2868 ... 


But 221930581 = 31 - 41 - 283-617, is not a prime. Again, these examples also fall into 
the category of isolated exotic expressions. The approximations are very precise but there 
are too few examples. So, we are limited to the few cases of primes appearing in these 2- 
term sums. But there is a trick, by taking the sum or the difference of 2 different values 
(different k). For example 


i nd n 
1 =-504) aq = 232256) say +6) oy =? aa 
n=1 


Who provides the approximation 


iets 61425 321300 


= —z— = 30.0299856... 


The approximations are less spectacular here, but what we gain is generality. By 
combining 4 terms of the type and with r = 1 or 4 and k odd. 


SOS. sae 
n=1 


one can easily obtain several representations for each prime even from 3. 2 having been 
obtained with Tangent numbers. Not only does it work quite well for all the primes up to 
1000000 but also one can choose an arbitrary prime. 


A test has been done with the first ones of the form 10” + c,,with m = 1 up to 101. These 
are the first prime numbers greater than 10” if m = 100 then c,, = 267. It is the sequence 
2,11, 101, 1009, ... (the sequence A003617 of the OEIS catalog). 


Example with m = 18, 


5778009767428887 783353934375 , 1246249957994342793434971153125 
= age a eg ee 


101843 = 
Which is 1000000676938336801.2703... 
The true value of the former is 


19 19 


(ee) A oo n 
1018 + 3 = 5937369506 >: om] — 6225783167123456 > eam =r 
n=1 n=1 
bad n23 oe n23 


n=1 n=1 
And for each of the terms of this sum we have the approximation (3), obviously the 
coefficients of the sums are prime to each other, two by two. 
k! 


CE aaa ~ S(k,r) 


which can be corrected with ¢(k + 1) to give an even more precise rational value if k is 
large, in other words 


k! 
S(k,r) Gayest) (8) 


which is rational when k is odd. More exactly when k = 4n +1. Whenk = 4n+3 it 
takes 2 terms in S(k,r). If we check the previous calculation we obtain the sum of 4 
rationals. 


67943139037162110976 518365013406083 


18 Ba ke 
eS 825 ss 6600 
248914919490662170624 118691882844287 
1365 10920 


There is something rather mysterious going on between equation (8) which gives a 
rational and the true value of S(k,r) which is irrational except ifr = 2 and k is of the form 
4n+3 (Ramanujan formula). Here what we found is that the formula (8) gives the true 
value if we take 2 terms in S(k,7r). Here is an example, 


We know that the following identity is true: 


1= 16) mq = -256) samy =a 
n=1 n=1 


But knowing that $(3,1) and S(3,4) are likely to be irrational taken separately. Our 
approximate formula (8) gives us S(3,1) ~ — and S(3,4) = — and if we add the 
coefficients (here 16 and -256), we get the true value of 1. Indeed : 


16 (=)- 256 ( re 1 
15 3840 


We can see that when k is small, here 3, the approximation is not very good and despite 
that, the equation becomes exact if we sum 2 terms. If k is larger, the approximation is 
much better and the result of the sum of 2 terms is again exact as with the case of p = 
1018 + 3 and k = 19 and 23. 


So, with a little program to find the coefficients of the sums S(k,1r) for each prime, we can 
easily reach 1000000 in a few minutes. As the size of the prime increases, there are more 
exact candidate sums and as many approximations. 


Given that the choice of the prime to represent seems arbitrary it is reasonable to 
conjecture that for any prime P we have the sum 


(oe) (oe) wk 
=0) a ee tcl aa +4) aay (7) 


n=1 n=1 
Here k, and k, are odd and different, a, b, c, and d are integers. Moreover, the first P is 
approximated by 


k,! k,! 
a: (rm)yet4 rv (r)k2+1 
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Where u and v are integers too. It remains to know what S(k,r) represents exactly. From 
the Ramanujan identity we can easily find those with r = 1 or 4. Bill Gosper had already 
explored some of these identities, including this one 
ed + << (301 + 210V27"/4 V7 + 90V2 73/4) 
> SES a0 + 71/4 4+ 120V7 + 90V2 7 ae 
n=1 @ 


This allows us to find the explicit values for r = 1, 2, 4 in S(k,r). 


k 


Skr=y, 4 Fay Value Approximation 
ae 11 we 1 6 
Cn 320 T(3/4)8 240 m* 
ae 1 7 1 3 
(32) 320 T(3/4)® 240 8 14 
fag 11 1 3 
Ge) 5120 T(3/4)® 240 128 14 
Maes are 120 
G1) 64 [(3/4)!2 ° 504 ra 
1 15 
ee 504 an 
3. ow 1 15 
ee) 22 T@/4y2 * 504 29 716 
a 363. 1 7! 
aD 2125 (3/4) 480 78 
- 3. ont 1 315 
(7,2) 2125 T(3/4)*6 480 24 18 
a 363. 1 315 
Oe 2175 T(3/4)** 480 212 776 
189 7° 1 9! 
sn) 28 [3/420 | 264 710 
1 2835 
mere 264 B10 
189 7 1 2835 
a) 28 (3/42 | 264 213 710 
See 393309 1° 691 11! 
oD 66560 T(3/4)24 65520 2 
ae 189° 691 155925 
G12) 66560 1'(3/4)24 65520 24 712 
ae 393309 1° 691 155925 
ae 222513 T(3/4)24 65520 216 712 
68607.” 1 13! 
se) 210 T(3/4)28 24 m4 


1 6081075 
eae 24 ros 
68607” 1 6081075 

S(13,4) 224 1(3/4)28 a 24 218 7714 
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Appendix 


Table of approximations of primes from 2. 


First Approximation Value 
196 
Z — 2.0121 
n* 
184275 2972025 
103 —__— + —___ 102.2651 
8 1° 478 
798525 2168775 
163 ee 160.9663 
8 1° 478 
5778009767428887783353934375 
1018 + 3 = ———s.99 1000000676938336801.2703 
1246249957994342793434971153125 
8 124 
Raw table 
Prime, S(k, 1), S(k, 4) [vector that makes the sum = 0 | Approximation 
Example 
7 = —228 S(6,1) + 58368 S(5,4) + 33 S(7,1) + 33792 S(7,4) 
SO 78929. PE ereO7 
84 78 8 110 
7, 7, 9, [-1, -228, 58368, 33, 33792], -4578525/ 4/Pi *8+95893875/ 8/ Pi ~10 
23, 5, 7, [-1, 63, 4032, 28, -7168], 61425/8/Pi ~61562275/4/ Pi *8 
29, 5, 7, [-1, 126, 8064, 24, -6144] , 61425/4/Pi ~64240975/ 2/Pi *8 
31, 5, 7, [-1, 504, 32256, -64, 16384], 61425/Pi “6- 321300/Pi *8 
31, 7, 9, [-1, -424, 108544, 66, 67584] - 4257225/ 2/ Pi *8+95893875/ 4/ Pi ~10 
41, 5, 7, [-1, 252, 16128, 16, -4096], 61425/2/Pi ~6+80325/Pi *8 
47, 5, 7, [-1, 315, 20160, 12, -3072], 307125/8/Pi ~6+240975/ 4/ Pi *8 
53, 5, 7, [-1, 378, 24192, 8, -2048], 184275/4/ Pi %6+80325/ 2/ Pi *8 
59, 5, 7, [-1, 441, 28224, 4, -1024], 429975/8/Pi %6+80325/4/ Pi *8 
71, 5, 11, [-1, 1890, 120960, -4, 16384], 921375/4/Pi ~6-638512875/ 4/ Pi ~12 
71, 5, 7, [-1, 567, 36288, -4, 1024], 552825/8/Pi “6- 80325/4/ Pi *8 
89, 5, 7, [-1, 756, 48384, -16, 4096], 184275/2/Pi ~6-80325/Pi *8 
89, 7, 9, [-1, -556, 142336, 99, 101376], -11165175/4/ Pi *84+287681625/ 8/ Pi ~10 
97, 3, 9, [-1, -498, 7968, 33, 33792], -11205/4/Pi 4+95893875/ 8/ Pi ~10 
97, 5, 7, [-1, 126, 8064, 152, - 38912], 61425/4/ Pi ~6+1526175/ 2/Pi “8 
103, 5, 7, [-1, 189, 12096, 148, - 37888], 184275/8/Pi ~6+2972025/ 4/ Pi *8 
109, 5, 7, [-1, 252, 16128, 144, - 36864], 61425/2/Pi %6+722925/Pi *8 
109, 7, 9, [-1, -36, 9216, 33, 33792], - 722925/4/Pi ~8+95893875/ 8/ Pi ~10 
113, 3, 9, [-1, -242, 3872, 33, 33792], -5445/4/Pi 4+95893875/ 8/Pi ~10 
113, 5, 7, [-1, 1008, 64512, -32, 8192], 122850/Pi ~6- 160650/ Pi *8 
127, 5, 11, [-1, 315, 20160, 2, -8192], 307125/8/Pi ~6+638512875/ 8/ Pi “12 
127, 5, 7, [-1, 441, 28224, 132, - 33792], 429975/8/Pi ~6+2650725/ 4/ Pi *8 
131, 3, 9, [-1, 46, -736, 33, 33792], 1035/4/Pi ~4+95893875/8/Pi ~10 
131, 5, 7, [-1, 126, 8064, 216, -55296], 61425/4/Pi ~642168775/ 2/Pi “8 
137, 3, 9, [-1, 142, -2272, 33, 33792], 3195/4/Pi “4+95893875/ 8/ Pi ~1O 
137, 5, 7, [-1, 189, 12096, 212, -54272], 184275/8/Pi *6+4257225/ 4/ Pi “8 
139, 5, 11, [-1, 3087, 197568, -6, 24576], 3009825/8/Pi “6- 1915538625/ 8/ Pi ~12 
139, 5, 7, [-1, 567, 36288, 124, -31744], 552825/8/Pi ~6+2490075/ 4/ Pi “8 
149, 3, 9, [-1, 334, -5344, 33, 33792], 7515/4/Pi “4+95893875/ 8/ Pi “LO 
149, 5, 7, [-1, 315, 20160, 204, -52224], 307125/8/Pi ~6+4096575/ 4/ Pi *8 
151, 3, 9, [-1, 366, -5856, 33, 33792], 8235/4/Pi “4+95893875/ 8/ Pi ~1O 
151, 5, 7, [-1, 693, 44352, 116, - 29696], 675675/8/Pi 6+2329425/ 4/ Pi *8 
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157, 3, 9, [-1, 462, - 7392, 33, 33792], 10395/ 4/ Pi “4+95893875/ 8/ Pi “10 

157, 5, 7, [-1, 756, 48384, 112, -28672], 184275/ 2/ Pi “6+562275/ Pi “8 

157, 7, 9, [-1, -428, 109568, 99, 101376], -8594775/ 4/ Pi “8+287681625/ 8/ Pi “10 
163, 3, 9, [-1, 558, -8928, 33, 33792], 12555/4/ Pi ~4+95893875/ 8/ Pi “10 

163, 5, 7, [-1, 819, 52416, 108, -27648], 798525/ 8/ Pi “6+2168775/ 4/ Pi “8 

167, 7, 9, [-1, -168, 43008, 66, 67584], -1686825/ 2/ Pi “8+95893875/ 4/ Pi “10 
167, 9, 11, [-1, 132, 135168, -8, 32768], 95893875/2/ Pi “10- 638512875/ 2/ Pi “12 
173, 3, 9, [-1, 718, -11488, 33, 33792], 16155/ 4/ Pi ~44+95893875/ 8/ Pi “10 

179, 3, 9, [-1, 814, -13024, 33, 33792], 18315/ 4/ Pi ~4+95893875/ 8/ Pi “10 

181, 3, 9, [-1, 846, -13536, 33, 33792], 19035/ 4/ Pi ~4+95893875/ 8/ Pi “10 

181, 5, 7, [-1, 1008, 64512, 96, -24576], 122850/ Pi ~6+481950/ Pi “8 

191, 5, 7, [-1, 756, 48384, 176, -45056], 184275/ 2/ Pi “6+883575/ Pi *8 
Bibliography 


[1] Encyclopedia of Integer Sequences, N.J.A. Sloane, Simon Plouffe, Academic Press , San 
Diego 1995. 

[2] The OEIS, Online Encyclopedia of Integer Sequences, sequences : sequences 
A000464, A000111, A000708, A006873, A000364. 

[3] Milton Abramowitz and Irene Stegun, Handbook of Mathematical Functions with 
Formulas, Graphs, and Mathematical Tables, New York: Dover Publications , 1972 
(original in 1964). 

[4] Richard J. Mathar, Table of Dirichlet L-Series and Prime Zeta Modulo Functions for 


Small Moduli http://arxiv.org/abs/1008.2547 


[5] Eisenstein series https://en.wikipedia.org/wiki/Eisenstein_series 


[6] Eric Weinstein, Eisenstein series (MathWorld): 
https: //mathworld.wolfram.com/EisensteinSeries.htm] 


[7] Apostol, Tom M. (7990). Modular Functions and Dirichlet Series in Number Theory 
(2nd ed.). New York, NY: Springer. ISBN 0-387-97127-0. 

[8] Bill Gosper, personal communication 

[9] Serre, Jean-Pierre, Cours d'arithmétique, Presses Universitaires de France, 1994, 


12 


